We have investigated the exclusive, radiative B meson decays to lf|(1430) 
in 89 x 10 6 BB events with the BABAR detector at the PEP-II storage ring. 
We measure the branching fractions B (B° -> if|(1430)°7) = (1.22 ± 0.25 ± 
0.10) x 10" 5 and B (J5+ -» if 2 *(1430)+7) = (1.45 ± 0.40 ± 0.15) x 10" 5 , 
where the first error is statistical and the second systematic. In addi- 
tion, we measure the CP-violating asymmetry A C p[B° -> iff (1430)°7] = 
-0.08 ±0.15 ±0.01. 
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We have investigated the exclusive, radiative B meson decays to 7^1(1430) in 89 x 10 6 BB events 
with the BABAR detector at the PEP-II storage ring. We measure the branching fractions B(B° — > 
7f 2 *(1430)°7) = (1.22± 0.25 ±0.10) x 10~ 5 and B(B+ -» if 2 *(1430)+7) = (1.45±0.40±0.15) x 10~ 5 , 
where the first error is statistical and the second systematic. In addition, we measure the CP- 
violating asymmetry Acp[B° -» if|(1430)°7] = -0.08 ± 0.15 ± 0.01. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



In the Standard Model (SM), flavor-changing neutral 
currents (FCNC) are forbidden at the tree level. For 
example, there is no direct coupling between the b quark 
and the s or d quarks. Effective FCNC are induced by 
loop (or "penguin" ) diagrams, where a quark emits and 
reabsorbs a W thus changing flavor twice. 

The discovery of B — > #* (892)7 decay [lj verified the 
existence of penguin processes. The same publication 
also reported evidence for B — > #2 (1430)7, later con- 
firmed by the BELLE collaboration |2(. Detailed knowl- 
edge about the decays to resonant modes with masses 
higher than #*(892), such as the B -> # 2 *( 1430)7 de- 
cay, will provide a better understanding of the inclusive 
b — > S7 branching fraction in terms of the sum over exclu- 
sive modes 3]. This is important because the compar- 
isons between the inclusive theoretical and experimen- 
tal rates place strong constraints on physics beyond the 
SM [4J. The measurement of the CP asymmetry, defined 

as A C p = r|j?Z^j + r(gI/fj > places a further stringent 
test on the SM, because the theoretical uncertainty in 
the non-perturbative hadronic effects cancels []|. 

This study is based on 81 fb _1 of data collected at the 
T(AS) resonance ( "on-resonance" ) with the BABAR detec- 
tor at the PEP-II asymmetric e+ (3.1 GeV)-e~ (9.0 GeV) 
storage ring, corresponding to 89 x 10 6 BB pairs. We 
have also collected a data sample of 10fb _1 at 40MeV 
below the T(4S) energy ("off-resonance"). 

The BABAR detector is described in detail elsewhere . 
Charged particle trajectories are measured by a five-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber. Photons and electrons are measured in the 
barrel and forward end-cap electromagnetic calorimeters, 
consisting of 6580 thallium-doped Csl crystals. 

Charged particle identification is provided by the en- 
ergy loss (dE/dx) in the tracking devices and by a ring- 
imaging Cherenkov detector (DIRC). A K/tt Cherenkov 
angle separation better than 4 standard deviations is 
achieved for charged tracks with momenta below 3 GeV/c. 

We use Monte Carlo (MC) simulations of events in 
the BABAR detector based on GEANT4 Q to optimize 
our selection criteria and to determine signal efficiencies. 
These simulations take into account variations of the de- 
tector conditions and beam backgrounds over the data- 
taking period. 

The #2 (1430) is reconstructed from three modes 
#1(1430)° ™> K+7r ~ and # 2 *(1430)+ -> #+7r°,#°7r + . 
#° mesons are reconstructed from the decay Kg — > 



7r + 7r~. Here and throughout this paper the charge- 
conjugate decays are included implicitly unless otherwise 
stated. 

A photon candidate is defined as a localized energy 
deposition well contained within the calorimeter accep- 
tance, —0.77 < cos6> < 0.96, where 9 is the polar an- 
gle with respect to the detector axis. It must have a 
lateral energy profile consistent with a photon shower 
and must be separated by 25 cm from all other show- 
ers, both neutral and charged. To suppress photons 
from Tr°(r]) decays, we veto any photon candidate that 
combines with another photon of energy greater than 
50 (250) MeV to form a 77 invariant mass in the range 
115 < M 77 < 155 (508 < M 77 < 588) MeV/c 2 . 

The 7T° candidates are reconstructed from pairs of pho- 
tons that have an energy above 50 MeV and an opening 
angle less than 36 degrees; the invariant mass of the two 
photons is required to be in the range 115 < M 77 < 
150 MeV/ c 2 . The candidate's momentum is recalculated 
with a 7T° mass constraint to improve the energy resolu- 
tion. 

The # ± and 7r ± track candidates are required to be 
consistent with originating from the e + e~ interaction 
point (IP); this requirement rejects tracks from beam- 
material and beam-gas interactions. A track is identified 
as a kaon if it passes through the DIRC radiators, and 
the detected Cherenkov photons are consistent in time 
and angle with a kaon of the measured track momentum. 
A charged pion is defined as a track that is not identified 
as a kaon or an electron, based on dE/dx and the ratio 
of the track momentum to the associated shower energy 
in the Csl calorimeter. 

The Kg candidates are reconstructed from two oppo- 
sitely charged tracks, having an invariant mass satisfying 
489 < M„ +7r - < 507MeV/c 2 . We require that the K° s 
candidate form a vertex that is displaced from the IP by 
at least 0.2 cm and lie in a direction from the IP consis- 
tent with the K s momentum. 

The #2(1430) candidate is required to have a Ktt in- 
variant mass within 120 (110) MeV/c 2 of the known Kf 
(#2 + ) mass 0. For the #+7r~ mode, we require that 
the two tracks are consistent with originating from a com- 
mon vertex. 

The B candidates are reconstructed by combining one 
#2*(1430) and one 7 candidate. To isolate the B meson 
signal, we use two kinematic variables. The first, AE, is 
defined as the difference between the reconstructed en- 
ergy of the B candidate and the beam energy, which is 
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known to high precision. The second is the beam en- 
ergy substituted mass (toes), which is defined as rn^g = 
V E b ca m -Ps, where E heam = y/s/2, p B = Pk*+P 1 with 
Pk* and p 7 representing the momenta of the if 2 and the 
photon candidates. For signal events, AE and mjgg" peak 
at zero and at the B meson mass, tob, respectively. For 
the modes containing a single photon candidate, namely 
K + ir~ and Kgir + , we adopt a technique that rescales 
the measured photon energy in the center-of-mass (CM) 
frame (denoted by asteroids) E* with a factor «, deter- 
mined for each event, such that E* K , + kE* = E heam 
in the rest frame of the T(AS); this improves the origi- 
nal TO ES ( m Es) resolution from 3.0 to 2.7MeV/c 2 . We 
retain B candidates with the invariant mass closest to 
the if 2 (1430) mass if we find multiple candidates with 
| AS | < 0.3 GeV and toes > 5.2 GeV/c 2 in the same event, 
which occurs in 3.1, 6.3, and 4.9% of the events for the 
K + ir~ , KgTT + and K + ir° modes, respectively. 

The background has two components, one of which 
includes combinatorial background from B decays and 
continuum qq production, where q can be a u, d, s or 
c quark, with the high-energy photon originating from 
initial-state radiation (ISR) or from 7r° and i] decays. 
These backgrounds are non-peaking in toes and AE. 

The second background contribution is from other res- 
onant B — > X s j modes, predominantly B — > if* (1410)7, 
and non-resonant B — > Ktt^j decays. We label these the 
"peaking" background, since these decays have Toes and 
AE distributions similar to the signal. In order to dis- 
tinguish the B — > if|( 1430)7 signal from the background 
decays, we examine the helicity- angle distributions. The 
helicity-angle Oh is defined as the angle of the K + or 
if ° in the rest frame of the if.* (1430) with respect to 
the flight direction of the if 2 (1430), measured in the B 
meson rest frame. These modes have different helicity- 
angle distributions: sin 2 Oh cos 2 Oh for if.* (1430), sin 2 Oh 
for if* (1410) and primarily sin 2 Oh for non-resonant de- 
cays assuming J = 1 for the spin of the if ir system. The 
non-resonant decays may have higher angular momentum 
contributions but the lowest possible angular momentum 
state is dominant; therefore, the helicity-angle distribu- 
tion for the non-resonant decay is assumed to be the same 
as that of the B — + if* (1410)7 decay. The systematic 
uncertainty associated with this modeling is studied and 
included in the measured branching fraction uncertainty. 

We exploit the difference in the event topology be- 
tween signal and continuum background to reduce the 
continuum contribution. To remove radiative Bhabha 
and e + e~ — > t + t~ events, we require that the ratio of 
second-to-zeroth order Fox- Wolfram moments of the 
event be less than 0.9. The distribution of the thrust 
angle (9t, defined as the angle between the direction of 
the photon candidate and the thrust axis of the rest of 
the event in the CM frame, is shown in Fig. Ufa). The 
rest of the event includes all the particles not used in the 
reconstruction of the B candidate. 



We train a neural network [9j with a combination of 
the thrust angle, the angle of the B meson candidate's 
direction with respect to the beam axis, the scalar sum 
of CM momentum of the rest of the event H3 (binned 
with 10° intervals ranging from parallel to anti-parallel 
relative to the photon momentum), sphericity, and the 
ratio of second-to-zeroth order Fox- Wolfram moments in 
the photon recoil system, which suppresses ISR back- 
ground. The neural network improves background sup- 
pression significantly. The distribution of the neural net- 
work output (NNO) is shown in Fig. ^b) for MC signal, 
MC continuum background and off-resonance data. 




|cos 6t| NNO 

FIG. 1: The cosine of thrust angle (a) and neural network 
output (b) distributions of the B — ► if* (1430)7 MC simu- 
lation (filled circles), the off-resonance data (line), and the 
continuum background MC (open circles). The vertical line 
indicates the cut value. 

The cuts on thrust angle and neural network output 
have been optimized for the best statistical significance; 
an iterative method of optimization is used to minimize 
correlations. The optimized cuts are | cos$t| < 0.95 and 
NNO> 0.55, as indicated in Fig. HI 

The signal yields are extracted using a simulta- 
neous maximum-likelihood fit of the Toes, AE and 
|cos0#| distributions. The fit is performed indepen- 
dently for each of the decay modes considered here. 

The signal toes and AE distributions are well de- 
scribed by an asymmetric resolution function ("Crystal- 
Ball" function [ll|), having an approximately Gaussian 
core plus a long tail due to the energy leakage from 
the calorimeter for the photon candidates. The peaking 
background is assumed to have the same toes and AE 
distributions as the signal. The continuum background 
is parameterized empirically by an ARGUS function 0] 
for toes and a linear function for AE. 

The cos Oh distribution of the signal has been parame- 
terized with siv?0hcos 2 0h — X(cos 4 0h — cos 6 0h), where A 
is a parameter determined from the Monte Carlo sample 
to account for the effect of the detector acceptance and ef- 
ficiency. The I cos Oh | distribution of the "non-peaking" 
background is parameterized by a linear combination of 
exponential and constant components. 
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TABLE I: The efficiency, fitted signal yield, significance, and 
measured branching fraction B(B — > ifj (1430)7) for each 
ifS(1430) decay mode. 



Mode 


e (%) 


Signal 


Significance (<r) 


B(10- a ) 


K+tt- 


6.4 


69 ± 14 


5.7 


1.22 ±0.25 ±0.10 


K° s n+ 


1.9 


29 ± 10 


.3.1 


1.69 ±0.59 ±0.16 


K + TT° 


1.9 


20 ±9 


2.2 


1.23 ±0.55 ±0.15 



Figures [21 and show the toes, AE 1 , and 
| cos 9h I distributions for the three modes in data; also 
shown are the | cos 6h | distributions of the candidates 
in the signal region, -0.15 < AE < 0.10 GeV and 
5.272 < toes < 5.288 GeV/c 2 . The signal as well as back- 
ground yields are allowed to vary in the fit. All the non- 
peaking background parameters are determined by the 
fit. The signal and peaking-background helicity-angle, 
Crystal-Ball width, and shape parameters are fixed to the 
MC expectations. The means of the signal toes and AE 
functions are fixed to the MC expectations, calibrated us- 
ing B — ► if* (892)7 candidates from MC simulation and 
data, while the peaking-background means are allowed 
to float due to their complex composition. The signal 
yields are given in Table [i] The signal significance has 
been evaluated from the change in the likelihood when 
the fit is repeated with the signal yield set to zero, in- 
cluding the systematic uncertainties, which are assumed 
to be normally distributed. The branching fractions are 
calculated with the assumption that the T(4S) decays 
equally to neutral and charged B meson pairs. 
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FIG. 2: Distributions of (a) ities, (b) AE, and (c) | cos 9h\ for 
B° Al(1430)°7, A" 2 *(1430)° -> K+tt~ candidates in data, 
and (d) j cosfn | in the signal region. The solid line shows the 
result of the fit to the data. The peaking (dashed-dotted line) 
and non-peaking (dotted line) background contributions are 
also shown. 

Figure 0] shows the Kir invariant mass distribution 
where the cut on this quantity has been relaxed. The 
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FIG. 3: Distributions of (a) mES, (b) AE, and (c) | cosOh | for 
the B+ -> iS'|(1430) + 7, A" 2 *(1430) + -h. K%-k+ candidates in 
data, and (d) Icos^hI in the signal region. The solid line 
shows the result of the fit to the data. The peaking (dashed- 
dotted line) and non-peaking (dotted line) background con- 
tributions are also shown. The corresponding distributions 
for 7fJ(1430) + — * K + tt° candidates are shown in (e-h). 



non-peaking background, estimated from the data out- 
side the signal region, has been subtracted. The invari- 
ant mass is fit with a relativistic Breit-Wigner function 
plus a first-order polynomial. There is a clear enhance- 
ment around 1.4GeV/c 2 in both the neutral and charged 
modes. 

We use the kaon charge to tag the flavour and measure 
the direct CP asymmetry in the decay rate as Acp[B° — > 
iq(1430)°7] =-0.08 ±0.15. 

The systematic error on the branching fraction for each 
mode is shown in Table ITT1 The total systematic uncer- 
tainty is computed as the sum in quadrature of the com- 
ponents. The yields returned from the fit are divided 
by the number of BB events and corrected for the ef- 
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FIG. 4: (a) 7f + 7r" and (b) K%n + and if + 7r° invariant mass 
distributions for the signal region (see text) after background 
subtraction. 



ficiency to obtain the branching fraction; therefore, the 
1.1% uncertainty on the S-counting measurement is in- 
cluded. The AE resolution is dominated by the photon 
energy resolution, which is determined from data using 
7T° and r\ meson decays with symmetric daughter pho- 
ton energies. The deviation in the reconstructed r\ mass 
from the nominal r) mass provides an estimate of the 
uncertainty in the measured single photon energy. The 
photon isolation and 7r° jr\ veto efficiency depend on the 
event multiplicity, and the effect is estimated by "embed- 
ding" MC-generated photons into both an exclusively re- 
constructed B meson data sample and a generic B meson 
MC sample. The photon and tt° efficiency uncertainties 
are determined from a comparison of the efficiencies in 
data and MC for e + e~ — * t + t~ events. The uncertainty 
in tracking efficiency is estimated from a sample of tracks 
well measured in the SVT. We estimate the uncertainties 
in the K s efficiency by comparing the data and MC dis- 
tributions of the momentum and flight distance. The 
efficiency for kaon and pion identification in the DIRC is 
derived from a sample of the decays D* + — > D°ir + , with 
D° -> K-TT+. 

Because the variables used for the neural network 
training are mostly calculated through the information 
from the rest of the event, we use 3155 fully recon- 
structed B — > Dir~ candidates in data, as well as sim- 
ulated B — ► Dtt~ events, as control samples. The pion 
in the B — ► Dir~ decay is treated like the photon in the 
B — » if| (1430)7 decay for the calculations of the event 
variables; the difference in the efficiency of the selection 
on the neural network output between data and MC is 
used as the systematic uncertainty. The systematic er- 
ror also includes the uncertainty in the invariant mass 
and width of iQ(1430) an d its sub-mode branching frac- 
tions jl|. 

We estimate the systematic error due to the fitting 
procedures as follows. For the shape parameters of 
toes, AE, and Icos^l distributions, we vary the pa- 
rameters in the fit within their errors from the MC ex- 
pectations. We also test the validity of the peaking- 
background |cos0y| probability density function (PDF) 
by mixing up to 20% J = 2 components and generating 



TABLE II: Fractional systematic uncertainties (%) in the 
measurement of B(B — > 1430)7). 



Uncertainty 


K+n~ 


K°n+ 


K+n° 


BB events counting 


1.1 


1.1 


1.1 


Photon and 7r° detection efficiency 


2.5 


2.5 


7.5 


Photon energy scale 


1.0 


1.0 


1.0 


Photon energy resolution 


2.5 


2.5 


2.5 


Photon isolation 


2.0 


2.0 


2.0 


7T°/?7 veto 


1.0 


1.0 


1.0 


Tracking efficiency 


1.6 


0.8 


0.8 


Kaon identification efficiency 


1.0 


- 


1.0 


Pion identification efficiency 


0.6 


0.6 




Kg efficiency 




.3.0 




Sub-mode branching fraction 


2.4 


2.4 


2.4 


K*(1A30) mass/width 


1.6 


1.0 


1.1 


Signal PDF parameters 


3.9 


5.8 


6.3 


Background modeling 


2.6 


2.9 


2.9 


Peaking-background modeling 


3.5 


4.9 


4.8 


MC statistics 


2.5 


3.2 


3.2 


Total 


8.4 


10.2 


12.6 



MC samples with different PDF parameterizations. We 
use the largest deviation in these tests as the systematic 
error of the signal yield. There is also a systematic error 
associated with the limited statistics of the signal MC 
sample. 

The particle-antiparticle asymmetry in the detector re- 
sponse, which includes 0.35% uncertainty for the tracking 
efficiency and 1.0% uncertainty for the charged particle 
identification, predominantly contributes to the system- 
atic uncertainty of the Acp measurement. The uncer- 
tainty in the estimate of nuclear interaction asymmetry, 
which arises from the different interaction probabilities 
of K+ and R- and of tt+ and 7r _ , is 0.20% 0. The 
total Acp systematic uncertainty is 1.1%. 

We have presented a measurement of the branching 
fraction for B° -> if 2 *(1430)°7 of (1.22 ± 0.25 ± 0.10) x 
10 -5 , which has a 5.7a statistical significance; this is in 
agreement with, but more precise than, previous exper- 
imental results. We observe a signal with a statistical 
significance of 3.8cx for B + — > K% (1430) + 7 and measure 
the branching fraction to be (1.45 ± 0.40 ± 0.15) x 10~ 5 , 
by combining the results from K s tt + and K + ir° modes. 
Both results agree with the theoretical predictions based 
on a relativistic form- factor model . The Acp is mea- 
sured to be— 0.08±0.15±0.01, thus no evidence of direct 
CP violation is observed. 
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